Effect of three-body elastic scattering on heavy quark momentum degradation in the 

quark-gluon plasma 
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Heavy quark drag coefficients in the quark-gluon plasma are evaluated in the perturbative QCD. It 
is found that for charm quarks the contribution from three-body elastic scattering is comparable to 
those from two-body elastic and radiative scatterings while for bottom quarks it becomes dominant. 
Using a schematic expanding fireball model, effects of both two-body and three-body scatterings on 
the transverse momentum spectra of heavy quarks produced in Au-I-Au collisions at center of mass 
energy ^/Snn = 200 GeV are studied. Results on electrons from resulting heavy meson decays are 
compared with available experimental data. 

PACS numbers: 12.38.Mh;24.85.-|-p;25.75.-q 



One of the most interesting observations in central 
heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) is the suppressed production of hadrons 
with large transverse momentum (pt)0,I2|- This phe- 
nomenon has been attributed to the radiative energy loss 
of partonic jets produced from initial hard scattering of 
incoming nucleons as they pass through the created dense 
partonic matter |^^|^. The same mechanism fails, how- 
ever, to explain a similarly large suppression of high pT 
charmed mesons observed through their decay electrons 
as a result of the dead cone effect associated with massive 
quarks 0,0- Furthermore, experimental data have in- 
dicated that charm quarks develop a substantial elliptic 
flow in non-central heavy ion collisions at RHIC 8] , con- 
sistent with the prediction of the quark coalescence model 
that assumes a thermally equilibrated charm quark dis- 
tribution in the quark-gluon plasma (QGP) 0. In both 
the Fokker-Planck approach 10] and the transport model 
[III IT^ . to reproduce the observed strong suppression 
of charm production at high pr and large charm ellip- 
tic flow requires a much larger charm quark two-body 
elastic scattering cross section than that given by the 
perturbative QCD. Such a large cross section would re- 
sult if colorless resonances are formed in charm quark 
scattering with light quark ,13.] . On the other hand, it 
was recently realized [lj| that two-body elastic scattering 
with the pQCD cross section causes a similar energy loss 
for the charm quark as the two-body radiative scattering 
|15|. and a large fraction of the observed suppression of 
charmed meson production at high px can be accounted 
for when both effects are included. 

Since the density of the partonic matter formed in 
heavy ion collisions at RHIC is large, ranging from about 
1 fm~'^ near hadronization to more than 10 fm~^ during 
the initial stage, three-body elastic scattering may also 
contribute to charm quark energy loss in the QGP. Pre- 
vious studies have shown that gluon 16] and quark [T^ 
three-body elastic scatterings are more efficient than two- 
body elastic scattering for the thermalization of initially 
produced partons. In the present paper, the effect of 
three-body elastic scattering on heavy quark {Q = c, b) 
momentum degradation in the QGP is studied. 
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FIG. 1: Topologically different diagrams for heavy quark 
(thick line) three-body elastic scatterings by quarks and anti- 
quarks (thin line) with different flavors. Wiggle lines denote 
the gluon. 



For the processes Qqq — > Qqq, Qqq Qqq, and 
Qqq Qqq with different light quark (q) and antiquark 
(q) flavors, there are three topologically different dia- 
grams in the lowest-order QCD as shown in Fig. While 
diagram (c) corresponds to one diagram, two diagrams 
are generated from diagram (b), corresponding to either 
the left or the right gluon (wiggle line) is flrst exchanged, 
and four diagrams are generated from diagram (a) by ex- 
changing the two gluons in all possible ways but keeping 
only one attached to the heavy quark. 

To ensure true three-body scattering, the intermediate 
quark in both diagrams (a) and (b) of Fig. n]must be off- 
shell, and this can be achieved by including its collisional 
width in the QGP and keeping only the real part of its 
propagator in evaluating these diagrams as in the treat- 
ment of three-body scattering in hadronic matter |l8j . 
The quark width in QGP is given by T/h = J2i{\-^i\'^)i 
where the sum is over all scattering processes with ]A^i]^ 
being their squared amplitudes after averaging over the 
spins and colors of initial partons and summing over those 
of final partons. The symbol (• • •) denotes average over 
the thermal distributions of scattered quarks and anti- 
quarks in the QGP and integration over the momenta of 
all flnal-state partons. 

We first consider the quark width due to two-body elas- 
tic and radiative scatterings as well as the inverse process 
of the latter. For heavy quark two-body elastic scatter- 
ing, there is only one t-channel gluon-exchange diagram 
for the process Qq — s- Qq or Qq Qq, while there are 
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in addition one s-channel heavy quark pole and one u- 
channel heavy quark exchange diagram for the process 
Qg — !■ Qg. Diagrams for heavy quark two-body radiative 
scattering are then obtained by adding an external gluon 
to above diagrams, leading to 5 diagrams for the process 
Qq Qqg or Qq Qqg and 16 diagrams for the process 
Qg — > Qgg. For the light quark, there are additional dia- 
grams besides those similar to the ones for heavy quarks, 
and all these diagrams are included in our calculations. 
In evaluating these diagrams, we remove the collinear 
singularity in t-channel diagrams by using the screen- 
ing mass m£) = gT |19| for the exchanged gluon, where 
g is the QCD coupling constant and T is the temper- 
ature of the QGP, and the infrared singularity in two- 
body radiative scattering by including the thermal mass 
nig = niu I \pl 19] for the radiated gluon. 

Using the QCD coupling cng = g'^ /Air = 0.3, appro- 
priate for the energy scales considered here, and also in- 
cluding the thermal mass for time-like gluons as well as 
the thermal mass ~ mu/x/G jl^ for time-like light 
quarks, we find that both heavy and light quark coUi- 
sional widths in the QGP increase almost linearly with 
quark momentum as well as the temperature of the QGP. 
For a momentum of 4 GeV/c, the widths are about 75 
MeV for charm quark and 49 MeV for bottom quark 
at T = 175 MeV and increase, respectively, to about 
130 MeV and 93 MeV at T = 350 MeV. These widths 
are mainly due to two-body elastic scattering with only 
about 25% from two-body radiative scattering and its 
inverse process. For the light quark, its width is about 
50% larger than that of charm quark. Using above cal- 
culated quark widths in the quark propagator, we find 
that three-body elastic scattering by quarks and anti- 
quarks with different flavors increases quark widths by 
at most 10% and can thus be neglected. We note that 
since the heavy quark width is smaller than the gluon 
thermal mass (Fq < Wg), i.e., the time between heavy 
quark collisions is longer than the time for radiating a 
thermal gluon, the destructive Landau-Pomeranchuck- 
Migdal (LPM) interference effect which is neglected 
in present study, is expected to be small for heavy quark 
two-body radiative scattering as shown in Ref . |2(]| . 

In the Fokker-Planck approach, the momentum degra- 
dation of a heavy quark in the QGP depends on its 
drag coefficient, which is given by averages similar to 
that for the quark coUisional width, i.e., 7(|p|,r) — 
(|M|2) - (|M|2p . p')/|p|2 [llij]. In the above, p and 
p' are, respectively, the momenta of the heavy quark be- 
fore and after a collision. In Fig. El we show the charm 
(upper panels) and bottom (lower panels) quark drag co- 
efficients as functions of their momentum in QGP at tem- 
peratures T = 200 MeV (left panels) and T = 300 MeV 
(right panels). It is seen that for charm quarks the con- 
tribution from three-body scattering by light quarks and 
antiquarks with different flavors (dash-dotted line) is gen- 
erally smaller than those from two-body elastic (dotted 
line) and radiative scatterings (dashed line). For bot- 
tom quarks, the most important contribution is from 
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FIG. 2: Charm (upper panels) and bottom (lower panels) 
quark drag coefficients as functions of their momentum in a 
QGP of temperature T = 200 MeV (left panels) or T = 300 
MeV (right panels). 



two-body elastic scattering while two-body radiative and 
three-body elastic scattering give smaller but comparable 
contributions. We note that the contribution from dia- 
gram (b) of Fig. n is larger than that from diagram (a) 
by about one order of magnitude and that from diagram 
(c) by about two orders of magnitude. 

If the two light quarks or antiquarks in Fig. ^ have 
same flavor, the number of diagrams is then doubled 
from interchanging final two light quarks or antiquarks. 
These exchange diagrams give same contribution as that 
due to the direct diagrams in Fig. ^ For the pro- 
cess Qqq — > Qqq with same quark and antiquark fla- 
vor, besides diagrams similar to those of Fig. the 
quark and antiquark can annihilate to a time-like vir- 
tual gluon, giving rise to additional five diagrams cor- 
responding to different ways a gluon is exchanged be- 
tween the heavy quark and the light quark, antiquark, 
and time-like gluon. The contribution from these annihi- 
lation diagrams is found to be small, and this allows us to 
neglect the interference between these diagrams and the 
diagrams in Fig. ^ Neglecting also the small interfer- 
ence terms between above direct and exchange diagrams 
for scattering by two identical quarks or antiquarks, the 
contribution from quarks and antiquarks with same fla- 
vor is then the same as that from quarks and antiquarks 
with different flavors. 

For heavy quark three-body elastic scattering involv- 
ing gluons, there are 36 diagrams for the process Qqg — > 
Qqg or Qqg Qqg, and 123 diagrams for the process 
Qgg Qgg- We have not been able to evaluate all these 
diagrams. Instead, we assume that these processes are 
also dominated by diagrams similar to diagram (b) of 
Fig. n This is also partially supported by the finding 
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that for the process qgg — > qgg involving a Hght quark 
the result obtained from including thermal and screening 
masses in the Parke formula |23 that takes into account 
many diagrams with topologies different from that of di- 
agram (b) of Fig. n indeed gives a much smaller con- 
tribution to the quark drag coefficient. Within this ap- 
proximation, the contribution from Qgg — > Qgg to heavy 
quark drag coefficients is slightly smaller while that from 
Qqg Qqg and Qqg Qqg is about a factor of three 
larger than that due to scattering by quarks and anti- 
quarks, essentially due to the associated color and flavor 
factors in these processes. Because of the importance of 
Qqg Qqg and Qqg Qqg, the heavy quark width 
is significantly increased by three-body elastic scatter- 
ing. Including this additional width in the intermediate 
heavy quark propagator in diagram (b) of Fig. pleads to 
a charm quark drag coefficient due to three-body elastic 
scattering (dash-dot-dotted line in Fig. ^ that is com- 
parable to that due to two-body elastic or radiative scat- 
tering, and a bottom quark drag coefficient that is now 
dominated by three-body elastic scattering. The result- 
ing total drag coefficients for heavy quarks due to both 
elastic two- and three-body scattering are shown by solid 
lines in Fig. [5] 

To see the effect of three-body elastic scattering 
on heavy quark momentum degradation in QGP, we 
consider central Au-|-Au collisions at center-of-mass 
energy ^/sNN — 200 GeV. The initial pT spec- 
tra of charm and bottom quarks at midrapidity are 
taken to be dNJdPpT = 19.2[1 -f {pT/&f]/{{l + 
^<t/3.7)12[1 + exp(0.9 - 2pt)]} and dNt/d^pr = 
0.0025 [l -I- {pT/lQf] exp(-pT/1.495), respectively, with 
Pt in unit of GeV/c. Both are obtained by multiply- 
ing the heavy quark pT spectra from p+p collisions at 
same energy by the number of binary collisions (~ 960) 
in Au-|-Au collisions. For bottom quarks, their pT spec- 
trum in p-Hp collisions is taken from the pQCD predic- 
tion of Ref.|23] as it is expected to be more reliable, while 
for charm quarks it is determined instead from fitting si- 
multaneously measured pr spectrum of charmed mesons 
from d-|-Au collisions |2J| and of electrons from heavy 
meson decays in p-l-p collisions. In obtaining the latter, 
heavy quarks are fragmented to hadrons via the Peterson 
frag mentation function D(z) = l/{z[l — l/ z — e/{l — z)]"^} 
|25| . where z is the fraction of heavy quark momentum 
carried by the formed meson, with e taken to be 0.02 for 
charm quarks and 0.002 for bottom quarks in order to 
reproduce the empirical fragmentation functions used in 
Ref.[2^. These heavy quarks are initially distributed in 
the transverse plane according to that of the binary colli- 
sion number, and their transverse momenta are directed 
isotropically in the transverse plane. 

For the dynamics of formed QGP, we assume that it 
evolves boost invariantly in the longitudinal direction but 
with an accelerated transverse expansion. Specifically, 
its volume expands in the proper time r according to 
V{t) = ttR{t)^t, where R{t) = Rq + a/2(T - tq)^ is the 
transverse radius with an initial value Ro—7 fm, to=0.6 



fm is the QGP formation time, and a = O.lc^/fm is the 
transverse acceleration ^26.] . With an initial temperature 
Ti = 350 MeV and using thermal masses for quarks and 
gluons, this model gives a total transverse energy compa- 
rable to that measured in experiments. The time depen- 
dence of the temperature is then obtained from entropy 
conservation, and the critical temperature Tc = 175 MeV 
is reached at proper time Tc ~ 5 fm. 

For a heavy quark with an initial transverse momen- 
tum po, time evolution of its mean transverse momen- 
tum (pt) can be obtained from the Fokker-Planck equa- 
tion, i.e., d{pT)/dt ~ —{'^{pt,T)pt)- Parametrizing 
the momentum dependence of heavy quark drag co- 
efficients shown in Fig. |3 by j{pT,T) w 7o(T)[l -t- 
apx], we then have d{pT)/dt « — 7o((pt) + o,{Pt)) ~ 
— 7o((pt) + cl{pt)'^) if we take {p^) « (pt)^, which is 
valid for high transverse momentum heavy quarks as 
considered here. The final mean transverse momentum 
of the heavy quark after passing through the expand- 
ing QGP is then given by {pr) — B/{\ — aB), where 
B = poexp(-/^^^ 7o(T)dT)/(l + apa) with r/ denoting 
the smaller of the time when the QGP phase ends and 
the time for the heavy quark to escape the expanding 
QGP, which depends on its initial transverse momentum 
and position. The final px spectra of heavy quarks are 
then obtained by averaging over their initial spatial and 
transverse momentum distributions. 

Besides fragmenting to mesons, heavy quarks produced 
in heavy ion collisions can also coalesce or recombine 
with thermal quarks in the QGP to form heavy mesons 

Il2l|. In both hadronization mechanisms, the momen- 
tum spectra of formed mesons are softer than those of 
heavy quarks. Instead of including both contributions, 
we use the fragmentation model in the present study for 
simplicity. Because of their smaller scattering cross sec- 
tions with hadrons 27] and lower hadronic matter den- 
sity, momentum degradation of heavy mesons in subse- 
quent hadronic matter is small 15] and is thus neglected. 
The heavy mesons produced from heavy quark fragmen- 
tation are therefore allowed to decay directly to electrons 
in order to compare with those measured in experiments. 

In the left panel of Fig. |3| we show our results for 
the initial (dash-dot-dotted lines) and final pr spectra 
of heavy quarks in Au-|-Au collisions at y^SNN = 200 
GeV. For both charm and bottom quarks, their final pr 
spectra become softer with the inclusion of more scat- 
tering processes: dotted lines for two-body elastic scat- 
tering only, dashed lines for both two-body elastic and 
radiative scatterings, dash-dotted lines for adding also 
three-body scattering by light quarks and antiquarks 
with different flavors, and solid lines for further including 
other three-body scattering processes. Although bottom 
quarks are negligible at low transverse momentum, they 
are important at high transverse momentum as a result of 
their smaller momentum degradation in QGP than charm 
quarks. Due to larger bottom quark mass, electrons from 
their decays become dominant at high px as shown in the 
middle panel of Fig. |21 The ratio of the electron px spec- 
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FIG. 3: Initial and final transverse momentum spectra of 
lieavy quarks (left panel) and their decay electrons (mid- 
dle panel) as well as the nuclear modification factor Raa 
for electrons (right panel) in central Au+Au collisions at 
= 200 GeV. 



trum from final heavy mesons to that from initially pro- 
duced ones, defined as the electron nuclear modification 
factor Raa, is shown in the right panel of Fig. 13 It is 
seen that the effect from two-body elastic scattering (dot- 
ted hne) is as important as that from two-body radiative 
scattering, similar to that of Ref.0]. The electron Raa 
including both contributions (dashed line) is, however, 
slightly above the experimental data from the PHENIX 
collaboration (filled triangles) and the preliminary 
data from the STAR collaboration (open squares) |29l |. 
Adding the contribution from charm quark three-body 
elastic scattering leads to an electron Raa that is in rea- 



sonable agreement with the measured one as shown by 
the solid line. If we include only charm quark three-body 
scattering by light quarks and antiquarks of different fla- 
vors, which is more reliably computed in present study, 
the resulting electron Raa is shown by the dash-dotted 
line which barely lies on the upper error bars of experi- 
mental data. 

Our results indicate that heavy quark three-body elas- 
tic scattering is important for understanding the elec- 
tron nuclear modification factor in heavy ion collisions 
at RHIC. The most important heavy quark three-body 
elastic scattering process involves scattering with a gluon 
and a light quark or antiquark in the QGP, and this 
contribution has been evaluated using the assumption 
that they are dominated by t-channcl gluon-exchange di- 
agrams similar to diagram (b) of Fig. ^ for three-body 
scattering by quarks and antiquarks with different fla- 
vors. Although we have shown that this is a valid ap- 
proximation for the latter process, its validity in heavy 
quark three-body elastic scattering involving gluons re- 
mains to be verified. If more accurate calculations indeed 
give as large a three-body contribution as shown here, 
then scattering involving more than three particles may 
also need to be considered. To evaluate the latter contri- 
bution poses, however, a major theoretical challenge. 
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